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ABSTRACT 


Energy  levels  below  2.2  Mev  excitation  in  114'116Cd 
14  6 

and  Nd  have  been  investigated  using  the  (n,n'y) 
reaction.  Gamma— ray  angular  distributions  were  measured, 
and  energy  levels,  spins  and  gamma-ray  multipole  mixing 
ratios  deduced.  Spin  assignments  and  gamma-ray  multipole 
mixing  ratios  were  determined  by  comparing  the  experi¬ 
mental  data  to  the  predictions  of  the  Statistical  Compound 

Nuclear  Model. 

114 

For  Cd,  levels  and  their  spins  have  been  confirmed 
at  1134 ( 0+) ,  1210 ( 2+) ,  1283(4+),  1365(2+),  1733(4+)  and 
1843  (2+)  keV.  New  spin  assignments  have  been  made  for  the 
levels  at  1861  (3+)  and  2049  (3  +  ).  The  state  at  1959  keV 
has  been  shown  to  be  consistent  with  a  spin  assignment  of 
3”. 

116 

For  Cd ,  levels  and  their  spins  have  been  confirmed 
at  1213  (2+)  ,  1220  (4  +  )  and  1381  (0+)  keV.  New  spin  assign¬ 
ments  have  been  made  for  the  levels  at  1644  (2  +  )  and 
1917  (3  +  )  keV.  The  state  at  1923  keV  has  been  shown  to 

be  consistent  with  a  spin-parity  assignment  of  3  . 

146 

For  Nd,  levels  and  their  spins  have  been  confirmed 
at  1043  (4  +  )  and  1189  (3-)  keV.  A  new  spin  assignment  has 
been  made  for  the  level  at  1376(1)  keV.  A  state  at 
1787  keV  has  been  proposed. 
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Two  types  of  large  sample  experiments: 
a)  The  pulsed  beam  time-of-f light  method  used 
by  Cranberg  and  Levin  (Cr56)  .  The  scatterer 
used  was  a  cylinder  of  approximately  one  mole 
of  56Fe.  b)  The  ring  shaped  scatterer  method 
as  used  by  Nishimura  zt  at.  (Ni65)  .  The 
scatterer  was  composed  of  346.6  gm  of  56Fe.  19 

a)  The  geometry  of  the  beam,  target  and  de¬ 
tector  for  a  small  sample  experiment.  Two 
different  types  of  y-ray  experiments  have  been 
done  using  this  geonetry.  b)  The  angular  dis¬ 
tribution  for  the  1220  513  keV  transition  in 

116Cd,  as  measured  in  this  work.  The  two 
curves  are  the  distorted  theoretical  distribu¬ 
tions  for  different  neutron  energies.  The 
En  =  00  curve  corresponds  to  a  nonaligned  ini¬ 
tial  state.  The  energy  dependence  of  the  theo¬ 
retical  distribution  is  a  function  of  the 
coefficients.  For  example,  the  En  =  °°  angular 
distribution  is  a  factor  of  1.25  higher  than  the 
En  =  1.7  MeV  distribution  at  0y  =  90°.  This 
factor  decreases  to  1.15  for  tne  En  =  2.8  MeV 
distribution.  c)  The  excitation  functions  for 
the  1220  keV  state.  The  experimental  data  is 
that  of  Gill  zt  at  (6:74)  which  was  taken  at 
0y  =  90°.  An  isotropic  y-ray  angular  distribu¬ 
tion  was  assumed  in  this  work.  The  correspond¬ 
ing  theoretical  excitation  function  is  denoted 
by  the  solid  curve.  The  dashed  curve,  which  in¬ 
cludes  angular  distribution  effects,  was  ob¬ 
tained  by  the  use  of  a  correction  factor,  as  de¬ 
scribed  above.  The  correction  of  excitation 
curves  due  to  angular  distribution  effects  is 
str  ight  forward  when  there  is  only  one  y-decay 
mode  available.  This  procedure  is  complicated 
when  there  is  more  than  one  y-decay  mode,  as 
one  needs  to  know  not  only  the  spins  and  pari¬ 
ties  of  the  states  involved,  but  the  branching 
and  multipole  mixing  ratios  of  the  y-decays. 

This  is  self  defeating,  as  this  is  usually  the 
information  that  is  being  deduced.  22 

Experimental  set-up  and  electronics.  28 
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4.  Typical  y-ray  spectra  for  the  Cd  isotopes, 

a)  114Cd.  The  numbers  above  each  peak  refer 
to  the  energies,  in  keV,  of  the  y-ray  transi¬ 
tions  between  the  Cd  levels.  Background 
peaks  are  denoted  by  B.  b)  116Cd.  The  inset 
shows  the  y-ray  decays  of  the  triplet  of  states 
near  1920  keV.  34 

5.  The  deduced  level  and  decay  scheme  of  114Cd. 

The  dotted  lines  refer  to  transitions  that 
were  too  weak  to  provide  useful  information. 
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31. 

Typical  y-ray  spectra  for  the  1 4 6Nd  experiment, 
taken  at  En  <_  2.67  MeV,  0  =  90°.  a)  The 
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scale.  b)  The  scatterer  "out"  spectrum.  The 
data  collection  time  for  this  spectrum  was 
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approximately  1/ 2  of  the  time  taken  for  the 
spectrum  shown  in  a) . 
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32.  The  deduced  level  and  decay  scheme  of  146Nd. 

The  dotted  lines  refer  to  transitions  that 
were  observed,  but  were  too  weak  to  provide 
useful  information.  The  spin-parity  assign¬ 
ments  are  the  result  of  the  present  and  pre¬ 
vious  works.  75 
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CHAPTER  1 


INTRODUCTION 

Several  years  ago  Scharff,  Goldhaber  and  Weneser  (Sc55) 
pointed  out  that  certain  regularities  exhibited  by  the 
first  and  second  excited  states  of  even-even  medium 
weight  nuclei  could  be  described  by  a  collective  vibra¬ 
tional  model  of  the  nucleus.  The  H4,116c^  isotopes 

is  in  the  transition  zone  of  vibrational  to  rotational 
nuclei.  If  these  nuclei  are  well  described  by  the  col¬ 
lective  vibrational  model,  their  low-lying  energy  level 
scheme  should  have  the  spin  sequence  0,2,  (024), 

3  ,  (0+2+3+4+6+) .  At  the  time  of  the  Scharff,  Goldhaber 

and  Weneser  paper,  however,  little  information  about  the 
levels  above  the  second  excited  state  was  known.  As  a 
result,  several  experimental  techniques  have  been  used 
in  the  past  few  years  to  study  these  nuclei. 

For  -*--*-4  /  H^cd f  the  results  from  (p,p?)  and  (a,a') 
experiments  (St68,  Ko69,  Ha63,  Co65,  De72)  established 
the  level  scheme  up  to  3.0  MeV.  The  one  and  two  phonon 
states  of  these  isotopes  were  identified  through  the 

measurement  of  the  charged  particle  angular  distributions. 

114 

The  y-ray  decay  scheme  of  the  cd  levels  was  first 

113 

established  through  the  Cd(n,y)  work  of  Groshev  et  al . 
(Gr62) .  Backlin  et  al .  (Ba66)  repeated  the  neutron  capture 


1 


2 


work  while  also  measuring  the  internal  conversion  spectrum. 

This  resulted  in  a  spin-parity  assignment  of  0+  to  the 

114  115 

Cd  level  at  1306  keV  excitation.  As  Cd  has  a  life¬ 
time  of  only  53  hours,  no  neutron  capture  work  has  been 

1“|  C 

done  to  investigate  the  y-ray  decay  scheme  of  Cd . 

Angular  distribution  measurements  of  y-rays  from 

the  114 ' 116Cd  (p ,  p  1  y )  [Mi69,  De73]  ,  114 ' 116Cd  (a , a ' y ) 

[Mi 6 9 ,  Ha63]  ,  and  114 ' 116Cd  (160 , 160 ' y )  [Mc65]  reactions 

corroborated  the  earlier  spin-parity  assignments  for  the 

one  and  two  phonon  states.  In  addition,  single  particle 

states  [1306  (0+)  and  1365  (2+)  keV  in  114Cd  and  1280  (2+) 
116 

keV  in  Cd]  were  observed  at  approximately  the  same 

energies  as  the  two  phonon  states. 

Investigations  of  the  excitation  functions  of  y-rays 

following  the  ^4  '  ^"^Cd  (n  ,n  '  y )  reactions  (Gi74)  revealed 

further  features  of  the  level  structures  of  these  two 

114 

nuclei.  The  states  at  approximately  1860  keV  m  cd 

116 

and  1920  keV  in  cd  were  found  to  be  a  doublet  and 

triplet  respectively.  Also,  the  previously  proposed 

114 

levels  at  1644  keV  excitation  m  Cd  and  1733  keV 

116  — 
excitation  in  cd  were  confirmed.  The  previous  3 

assignments  to  the  ^4  >  H^cd  states  near  1960  keV  and 

1920  keV  were,  however,  disputed.  Their  results, 

together  with  the  results  from  previous  works  are  listed 

in  Tables  I  and  II. 

The  earliest  works  resulting  in  spin-parity  assign- 

146 

ments  for  the  excited  states  of  Nd  involved  Coulomb 


I 


Table  I.  A  Summary  of  the  Results  of  Previous 

114 

Investigations  of  Cd. 


Level  Energy  (keV) 

(p,p')a) 

(n,y)b) 

,  ,  .d) 

(n,n'y) 

558 

2  + 

2  + 

2 

1134 

0  + 

0  + 

0 

1210 

2+ 

2  + 

2 

1283 

4  + 

4  + 

4 

1306 

0+c) 

0 

1365 

2  + 

2  + 

2 

1733 

4  + 

4 

1843 

2 

1 

1861 

(0) 

1865 

(0+) 

1 

1959 

3“ 

( 2~) 

(1) 

2049 

1+ 

(4) 

a) 

Ref . 

Ko69 

b) 

Ref . 

Gr68 

c) 

Ref . 

Ba65 

d)  Ref.  Gi74 
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Table  II.  A  Summary  of  the  Results  of  Previous 

116 

Investigations  of  Cd. 


Level  Energy  (keV) 

(n,  n  '  y) 

Excitation 

Previous 

Results 

513 

2 

2+d> 

1213 

2 

2b) 

1220 

4 

4b) 

1283 

0 

2C) 

1381 

0 

0d) 

1644 

(2,3) 

1917 

(0) 

1923 

3“  (~  1920 ) d) 

1930 

2120 

(3',5_)  (^2115)d) 

a) 

Ref . 

Gi74 

b) 

Ref . 

De7  3 

c) 

Ref . 

De7  2 

d) 

See 

Table  1  of  Ref.  De72  and  references  therein 

V 
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excitation  (Ha62,  Be67,  Bu67) .  Spin-parities  of  2+,  4+, 

and  3  were  assigned  to  the  levels  at  453,  1043  and  1189 

keV  excitation  on  the  basis  of  y-ray  angular  distribution 

measurements.  The  y-ray  decay  scheme  for  higher  levels 

was  established  through  (n,y)  work,  as  compiled  by 

1 4 

Groshev  at  at.  (Gr68) ,  The  level  scheme  of  Nd  was 
also  investigated  by  the  (d,d')  work  of  Christensen  at 
at.  (Ch70)  and  the  (t,p)  work  of  Chapman  at  at.  (Ch72) . 

In  those  works,  spin-parity  assignments  made  on  the 
basis  of  the  angular  distributions  of  the  outgoing 
charged  particles  were  consistent  with  those  of  earlier 
works.  These  results  are  listed  in  Table  III. 

Although  many  experiments  have  been  carried  out  in 
an  effort  to  investigate  the  structures  of  these  three 
nuclei,  a  few  uncertainties  remain.  The  situation  re¬ 
garding  spin  assignments  for  the  levels  above  1300  keV 
excitation,  and  the  multipole  mixing  ratios  for  the 
y-ray  decays  is  still  uncertain. 

Neutron  scattering  reactions  represent  a  valuable 
technique  for  such  investigations  since  the  absence  of 
a  Coulomb  barrier  allows  the  levels  of  interest  to  be 
excited  close  to  threshold,  thereby  eliminating  compli¬ 
cations  due  to  cascade  feeding.  In  this  work  y-ray  ang¬ 
ular  distributions  from  the  '  ^^Cd  (n ,n 1  y )  reactions 
have  been  studied,  providing  further  information  about 
the  spins  and  parities  of  the  low-lying  levels. 
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Table  III.  A  Summary  of  the  Results  of  Previous 

.  14  6 

Investigations  of  Nd. 


Level  Energy  (keV) 

(160,160'y)a) 

(d,d')b) 

, ,  x  c) 

(t,p) 

454 

2  + 

2  + 

(2+) 

1043 

4  + 

4+ 

1189 

3" 

3“ 

1377 

(1“)  -  1373 

1471 

(2+)  ~  1468 

1515 

(5”) 

1680 

0+  ~  1693 

1740 

4  + 

a) 

Ref . 

Bu67 

b) 

Ref . 

Ch7  0 

c)  Ref.  Ch72 
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Multipole  mixing  ratios  for  the  reported  transitions 
were  also  deduced.  This  technique  has  the  following  ad¬ 
ditional  advantages: 

1.  The  three  phonon  states  near  1900  keV  can  be 
investigated,  which  was  not  possible  using 
Coulomb  excitation  (5-9) . 

2.  There  is  no  necessity  of  detecting  all  of  the 
y -decay  modes  from  a  given  level. 

3.  The  systematic  error  associated  with  y-ray  ex¬ 
citation  functions  due  to  the  angular  depend¬ 
ence  of  the  y-ray  is  avoided. 

4.  The  sensitivity  of  the  angular  distribution  cal¬ 
culations  to  the  optical  model  parameters  is 
much  smaller  than  it  is  for  excitation  functions. 


These  last  three  points  were  a  source  of  error  in  the 
(n,n'y)  excitation  work  of  Gill  dt  at. 


■ 
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CHAPTER  2 

THEORY 

In  the  compound  nucleus  picture  of  a  nuclear  reaction, 
the  energy  carried  into  the  target  by  the  projectile  is 
shared  with  all  of  the  nucleons  shortly  after  contact. 

Once  the  energy  is  shared,  a  large  number  of  energy  ex¬ 
changes  are  required  for  one,  or  a  system  of  particles  to 
acquire  enough  energy  so  that  it  can  be  emitted  by  the 
compound  nucleus.  Because  of  the  energy  sharing,  the  dis¬ 
integration  of  the  compound  system  does  not  depend  on  the 
method  by  which  it  has  been  formed.  The  mode  of  decay  does, 
however,  depend  on  the  constants  of  motion:  energy,  angular 
momentum,  and  parity.  As  a  result,  it  is  possible  to 
separate  a  nuclear  reaction  into  two  separate  events;  one, 
the  formation  of  the  compound  nucleus,  and  two,  the  dis¬ 
integration  of  the  system  into  the  products  of  the  reaction. 
This  is  the  basis  of  the  compound  nucleus  model. 

To  determine  the  cross  section  for  the  formation  of 
the  compound  nucleus,  all  possible  two-body  interactions 
between  the  projectile  and  the  nucleons  of  the  target  should 
be  considered.  However,  this  approach  is  much  too  compli¬ 
cated  except  for  the  simplest  reactions.  Consequently,  a 
model  in  which  the  many  body  problem  is  reduced  to  that  of 
two  structureless  bodies  acting  through  an  average  potential 


V 


is  adopted.  This  approximation  is  known  as  the  optical 
model . 

In  the  optical  model,  the  incident  particle  is  in¬ 
fluenced  by  a  potential  of  the  form 


u (r )  =  Vc(r)  -  v0f(r,R0,aQ)  -  i (w  -  4a;jWD|-J f (r ' Rk 


+2(vso  +  iwso)  Fa-  f(r,R__,a„)  3'  £ 


so  so 


( 


where  Vq  is  the  real  well  depth,  W  and  WD  are  the  volume 

and  surface  imaginary  well  depths,  and  V  and  W  are 

s  o  s  o 

the  real  and  imaginary  parts  of  the  spin-orbit  well  depths 

The  functions  f(r,R  ,a  )  have  the  form 

x  x 

IT“"R 

f(r'Rx,ax)  =  d  +  exp  )  _1  ( 

x 

where  R  and  a  are  the  corresponding  radii  and  well 

X  X 

diffuseness  parameters.  Vq (r)  is  the  Coulomb  potential 

which  is  taken  to  be  that  due  to  a  charge  of  Ze  uniformly 

distributed  over  a  sphere  of  radius  .  The  real  and 

c 

imaginary  potentials  account  for  refraction  and  absorp¬ 
tion  of  the  incident  particle  whereas  the  spin-orbit 
potential  term  takes  into  account  their  polarization. 

Once  formed,  the  compound  nucleus  can  decay  to  a 
final  state  through  what  is  called  a  reaction  channel. 

The  final  states  of  the  nuclear  reaction  can  be  linked 
to  the  initial  state  by  a  collision  matrix  U.  The  colli¬ 
sion  matrix  is  composed  of  elements  which  relate  the 
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amplitudes  of  the  various  partial  waves  of  the  reaction 
channels  to  the  partial  waves  of  the  incident  channel. 
This  can  easily  be  seen  if  one  considers  the  total  wave 
function  of  the  nuclear  reaction  ¥ .  It  can  be  written  as 


E  v-l/2 
c  c 


-  BO  ) 

C  V 


(3) 


where  the  summation  is  over  all  allowed  channels.  The 

parameters  I  and  O  represent  the  incoming  and  outgoing 
c  c 

waves,  having  amplitudes  of  Ac  and  Bc  respectively.  The 

factor  v  1/2  is  used  as  v  I  corresponds  to  the 
c  c  c 

incoming  flux.  \p  is  the  channel  wave  function,  which 
contains  the  quantum  description  of  the  channel.  For  a 
single  entrance  channel  c,  the  relationship  between  the 
amplitude  coefficients  is  given  by 


B 


c ' 


U 

c 


A 

c  c 


(4) 


where  U  .  is  an  element  of  the  collision  matrix.  U, 
c  c 

although  it  is  expressed  through  the  asymptotic  behavior 

of  the  wave  functions,  does,  however,  contain  the 

inherent  physics  of  the  collision  interaction.  By 

hypothesis,  each  matrix  element  can  be  decomposed  into 

a  weakly  energy  dependent  U  ,  and  a  rapidly  fluctuating 

o  c 

term  U  .  that  averages  to  zero  over  a  large  enough 
c  c 


energy  interval . 


That  is. 
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c 

o 

o 

— 

u  ,  +  u  , 

c  '  c  c  c 

where 

U  , 

c '  c 

= 

<U  ,  >  =  const 
c  c 

and 

*N» 

U  , 

c '  c 

7* 

const  but  <U  ,  >  =  0 
c  c 

Using  this  notation,  the  formation  of  the  compound 
nucleus  through  channel  c  can  be  written 


-  _  v2  I 

°c  -  kc  c 


<  I  u 


> 


c  c 


(1 


c 


T 

c 


(5) 


where  T  is  the  transmission  coefficient  for  channel  c, 
c 

and  k  is  the  wave  number  for  the  relative  motion  of  the 
c 

particles  in  channel  c.  The  transmission  coefficient 
expresses  the  flux  associated  with  the  compound  nucleus 
component  of  the  reaction.  The  total  cross  section  for 
a  reaction  channel  c*  is  the  product  of  the  factor  a 
that  expresses  the  cross-section  for  formation  through 
channel  c  and  a  branching  ratio  for  the  decay  of  the 
compound  system  into  exit  channel  c'.  Through  time 
reversal  invariance  and  conservation  of  flux,  the  problem 
of  the  decay  of  the  compound  nucleus  through  channel  c ' 
can  be  related  to  the  situation  where  the  compound 
nucleus  is  formed  through  channel  c'.  That  is, 


k2  5 


cc  1 


=  k2  a 

c 


c '  c 


(6) 


Use  of  the  principle  of  detailed  balance  allows  the 


average  cross  section  to  be  expressed  as  a  function  of 
the  transmission  coefficients. 


c 


c  1  c 


T  T  , 
IT  C  C 


k2  Z  T 

C  _ it 


where  the  summation  is  over  all  possible  exit  channels  c' 
The  labels  c,  c',  and  c'1  characterize  the  quantum 
description  of  the  respective  channels,  and  thus  contain 
a  prescription  of  the  various  angular  momenta  and 
parities  involved.  When  this  dependence  is  explicitly 
taken  into  account,  one  finds  (Vo73)  that  when  a  compound 
nucleus  is  formed  through  channel  c,  the  average  cross 
section  for  the  decay  through  channel  c'  is 

ac'c  “  k2  J7T  (21  +  1)  ( 2  i  +  1 )  1  j£  3 


x{. 

3 


Z 

V 


Z 

.  t  * 


z 

j '  ' ' 


where  I  and  i  are  the  respective  spins  of  the  target  and 
projectile.  This  is  the  basis  of  the  statistical  theory 
of  nuclear  reactions. 

The  exit  channels  under  consideration  in  this  work 
are  those  that  involve  gamma  radiation.  For  the  gamma 
decay  from  the  initial  state  | of  the  residual 
nucleus  to  the  final  state  the  probability 

amplitude  is  given  by  (Ro67) 


' 
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1M2 


(k) 


k  ,1/2 
2iTh  i 


I 

LMtt 


IT 

q 


<jimi 


1<TT>| 

'LM1 


J2M2> 


(R) 

Mq 


(9) 


where  T  represents  the  interaction  in  terms  of  the 

electromagnetic  operators;  <tt  =  0 >  denoting  the  electric 

operators  and  <tt  =  1>  the  magnetic  operators.  The 

quantity  q  takes  on  the  values  of  ±1  for  right-handed  and 

left-handed  circular  polarization  respectively.  The 

rotation  matrix  DL  (R)  describes  a  rotation  that  projects 

Mq 

the  propogation  vector  £  onto  the  z-axis.  When  the  spin 

orientation  of  the  final  state  is  not  observed,  the 

transition  probability  Pq  (k)  for  a  gamma  emission  from 

M1 

the  initial  state  is  an  incoherent  sum  over  the  possible 
values  of  .  The  total  transition  probability  for 
photons  with  the  momentum  k  is  then  obtained  by  weight¬ 
ing  each  Pq  (k)  with  the  relative  population  parameters 
w(M^)  of  the  substate  M  . 


pq  (k) 


(Mi) 

w  (M  ) 
1  1 


(k) 


Aq 

M]m2 


(10) 


The  angular  distribution  for  the  emission  of  photons  from 
the  state  to  the  state  J2  is  obtained  by  an  expansion 
of  (10) 


w(0)  =  pq(k)  =  27h  LL'K  ttu'  VA1  Rk(LL'J1J2>  PK(COS0>  x 

<TT  > 


TT  +  TT  1  <J. 

q  X  1 


!  t<7T>Ii  j^> 


<J. 


T 


L’ 


J2>* 


(ID 


( 2L+1 ) 


1/2 


( 2L 1 +1 ) 


1/2 


. 
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where 

BK(J1}  =  M  W(M1)  (-1) J1"M1  (2J1+1)1/2  ( J1 1 KO)  , 
and 

RgtLL'J^)  =  (-l)q+Jl"J2+L'-L-K  (2J1  +  1,1/2  (2L+1)1/2 

1/2 

x(2L'+l)  (LL  '  q-q  |  KO)  W(J ^  I*'  jKJ 2)  . 


The  P^(cos0)  are  the  Legendre  Polynomials.  As  the 
coefficients  are  functions  of  the  population  parameters, 
they  describe  the  effects  due  to  nuclear  alignment.  The 
R^(LL'J^J2)  coefficients  depend  only  on  the  quantities 
that  characterize  the  gamma  transition.  The  reduced 

l/2 

matrix  elements,  divided  by  factors  of  (2L+1)  ,  are 

proportional  to  the  square  root  of  the  partial  gamma 
width  due  to  the  multipole  transition  (tt,L).  That  is; 


4k 


ttL 


<J 


<TT  > 


1 


J2> 


/  ( 2L+1 ) 


(12) 


The  reduced  matrix  elements  contain  the  nuclear  informa¬ 
tion.  Experimentally,  the  sum  of  their  square  moduli 

£  r  \ 

can  be  obtained  by  lifetime  information  (1/t  =  (1/h)  ttL'  * 

ttL 

Consequently,  for  angular  distribution  measurements, 
one  is  interested  in  determining  only  the  ratios  of  the 
reduced  matrix  elements.  To  further  this  end,  the  multi¬ 
pole  mixing  ratios  are  introduced. 

6<*>  .  <JillT^>||J2>/(2L+l)V2 
L  <J]_ |  iTj/Z  |J2>/(2L+1)1/2 


(13) 


15 


where  L  tt  represents  the  lowest  possible  multipolarity 
involved  in  the  gamma  decay.  Using  the  mixing  ratio, 
and  summing  over  q,  as  polarization  was  not  observed  in 
this  work,  Eq.  (11)  becomes 


W(9)  =  Z  Z  {B  ( J-i  )  Rt.(LL'J.J0)  P^(cos0)} 
LL '  Ktttt  k  1  K  1  2  K 


x  {  l+(-l)L+L,+7I+,T,-K  }  6  <7T>  s™'  >  /  2  I  6  <1T>  I  2 

L  L'  LTf  L  1 


The  summation  is  over  only  K=even  terms  when  the  nucleus 
is  in  an  aligned  condition  [w(M^)  =  wC-M^)];  therefore 
Bt,  __  =  0).  This  same  condition  holds  true  even  when 
the  initial  state  is  not  aligned,  if  and  J2  have 
definite  parity.  As  electromagnetic  operators  conserve 
parity,  there  are  selection  rules  for  the  type  of  radia¬ 
tion  that  can  occur  for  a  gamma  transition  between  states 
of  parity  1^  and  n  .  They  are 


(n1)  <n2) 

Type  of  2 

tt=0  (electric) 

radiation 

tt=1  (magnetic) 

1 

even 

odd 

-1 

odd 

even 

Hence,  all  terms  in  Eq.  (14)  have  L+L  ' +tt+tt  '  =even ,  and 
the  factor  (1  +  (-1)L+L  +7T+Tr  K)  are  zero  if  K  is  odd, 
and  2  if  K  is  even.  If  only  two  multipoles  contribute 
to  the  transition,  the  theoretical  angular  distribution 


(14) 


becomes 
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K=even  BK(J1>  pK(cos(e)  x {VLLJ1J2>  + 
26Rk(Elj1J2)  +  62Sk(LLJ1J2)}/(1+62) 


This  equation  is  valid  when  the  two  states  have  definite 
spins  and  parities  and  when  circular  polarization  is  not 
observed.  An  experimental  angular  distribution  from  an 
aligned  nucleus  has  the  form 


W 

exp 


(0) 


£ 

K=even 


aK  PK(cos0) 


(16) 


To  facilitate  the  comparison  of  the  experimental  data  to 
the  theoretical  angular  distribution.  Green  et  al .  (Gr72) 

have  rewritten  Eq.  (15)  in  the  form  of  Eq.  (16)  with  the 
transformation 

aK  =  aK  +  eKsin2(p-pK0) 

where  p  =  arctan6  (17) 


The  coefficients  a^,  3^,  and  are  functions  of  the 

and  coefficients. 

In  principle,  if  the  B^  coefficients  are  known,  it 
is  possible  to  extract  the  multipole  mixing  ratio  6  from 
an  experimental  y-ray  angular  distribution.  The  mixing 
ratios  are  defined  in  terms  of  the  reduced  matrix  elements, 
which  in  turn  are  functions  of  the  matrix  elements 

<J1M1 ' TLM  lJ2M2>*  Therefore,  by  deducing  the  value  of  5, 
information  about  the  initial  and  final  state  nuclear 


wave  functions  can  be  obtained. 
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CHAPTER  3 

EXPERIMENTAL  PROCEDURE 

3 . 1  Data  Analysis 

The  (n,n'Y)  reaction  is  useful  in  nuclear  physics 
as  it  enables  one  to  investigate  the  structure  of  a  heavy 
target  nucleus  at  a  low  energy.  Reactions  initiated 
with  charged  particles  are  not  feasible  at  these  low 
energies  because  of  the  large  Coulomb  barrier.  Ideally, 
a  neutron  scattering  experiment  would  utilize  a  small 
sample  situated  at  a  large  distance  from  the  neutron 
production  source.  This  would  guarantee  that  all  reac¬ 
tions  in  the  scatter  would  be  initiated  by  monoenergetic 
neutrons  colinear  with  the  beam  axis.  However,  this  is 
not  feasible  because  of  the  difficulty  in  obtaining  a 
high  flux  neutron  beam.  To  reduce  the  time  involved  in 
obtaining  a  reasonable  amount  of  data,  the  large  sample 
method  (  Di72  and  references  therein)  was  formulated 
(see  fig.  1) .  It  required  a  large  amount  of  scatterer 
material  (approximately  one  mole)  being  placed  at  a  large 
distance  away  from  the  neutron  production  source.  This 
had  the  effect  of  maximizing  the  number  of  reactions 
while  restricting  the  energy  spread  of  the  neutron  beam 
to  an  acceptable  level.  For  most  samples  however,  the 
target  cost  is  prohibitively  high.  This  resulted  in  the 
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Figure  1 . 


Two  types  of  large  sample  experiments:  a)  The  pulsed  beam 
time-of-fl  ight  method  as  used  by  Cranberg  and  Levin  (056). 
The  scatterer  used  was  a  cylinder  of  approximately  one  mole 
of  56Fe.  b)  The  ring  shaped  scatterer  method  as  used  by 
Nishimura  et  al .  ( N i 6 5 ) .  The  scatterer  was  composed  of 
<46.6  gm  of  56Fe. 
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development  of  the  small  sample  method  (E171)  for  study¬ 
ing  nuclei.  The  small  sample  technique  relied  on  a  close 
scatterer-neutron  source  geometry.  Although  this  geometry 
gave  rise  to  a  high  neutron  flux,  there  was  a  large 
spread  in  neutron  energy,  corresponding  to  the  differing 
angles  with  which  the  neutrons  left  the  production  source. 

Another  inherent  difficulty  involved  with  the  close 
geometry  method  is  the  distorting  effect  of  the  finite 
sizes  of  the  scatterer  and  the  neutron  production  source. 
Corrections  to  y-ray  angular  distribution  data  to  account 
for  these  effects  can  only  be  made  if  the  spins  and 
parities  of  the  states  involved,  and  the  multipole  mix¬ 
ing  ratio  of  the  y  decay  are  known.  As  this  is  usually 
the  information  that  is  desired,  the  distorting  effects 
must  be  folded  into  the  theoretical  angular  distributions. 
The  distorted  theoretical  angular  distributions  can  then 
be  compared  to  the  experimental  data  to  deduce  the  y-ray 
multipole  mixing  ratio,  and  the  initial  state's  spin  and 
parity.  This  method  was  developed  by  Davidson  et  al . 

(Da7  6a) . 

The  geometry  of  the  experiment  is  shown  in  figure  2. 

The  incident  proton  of  energy  E  reacts  with  the  neutron 

production  source  at  a  point  P  on  a  circular  beamspot  of 

radius  r.  A  neutron  of  energy  En  is  produced,  leaving 

the  production  source  at  an  angle  0n  with  respect  to  the 

proton  beam  axis.  The  neutron  then  travels  a  distance 

£  inside  the  scatterer  to  reach  the  point  P*  where  it 
n 


I  ; 


Figure  2. 

a)  The  geometry  of  the  beam,  target  and  detector  for  a 
small  sample  experiment.  Two  different  types  of  y-ray  ex¬ 
periments  have  been  done  using  this  geometry. 

b)  The  angular  distribution  for  the  1220-*513  keV  transi¬ 
tion  in  1 1 ^Cd ,  as  measured  in  this  work.  The  two  curves 
are  the  distorted  theoretical  distributions  for  different 
neutron  energies.  The  En  =  00  curve  corresponds  to  a  non- 
aligned  initial  state.  The  energy  dependence  of  the  theo¬ 
retical  distribution  is  a  function  of  the  coefficients. 
For  example,  the  En  =  00  angular  distribution  is  a  factor 
of  1.25  higher  than  the  En  =  1.7  MeV  distribution  at 

6y  =  90°.  This  factor  decreases  to  1.15  for  the  En  =  2.8 
MeV  distribution. 

c)  The  excitation  functions  for  the  1220  keV  state.  The 
experimental  data  is  that  of  Gill  et  al  (6:74)  which  was 
taken  at  0y  =  90°.  An  isotropic  y-ray  angular  distribu¬ 
tion  was  assumed  in  this  work.  The  corresponding  theoret¬ 
ical  excitation  function  is  denoted  by  the  solid  curve. 

The  dashed  curve,  which  includes  angular  distribution 
effects,  was  obtained  by  the  use  of  a  correction  factor, 
as  described  above.  The  correction  of  excitation  curves 
due  to  angular  distribution  effects  is  straight  forward 
when  there  is  only  one  y-decay  mode  available.  This  pro¬ 
cedure  is  complicated  when  there  is  more  than  one  y-decay 
mode,  as  one  needs  to  know  not  only  the  spins  and  parities 
of  the  states  involved,  but  the  branching  and  multipole 
mixing  ratios  of  the  y-decays.  This  is  self  defeating,  as 
this  is  usually  the  information  that  is  being  deduced. 


Beam  Production  Target 
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inelastically  scatters  from  a  target  nucleus.  The 
resulting  y-ray  is  emitted  in  a  direction  0^  with  respect 
to  the  proton  beam  and  if with  respect  to  the  incident 
neutron  direction.  If  the  y-ray  is  not  scattered  or 
absorbed,  it  travels  a  distance  of  Z^  before  leaving  the 
scatterer.  A  detector  placed  in  the  direction  of  0 
with  respect  to  the  beam  line  detects  y-rays  from  all 
possible  combinations  of  the  points  P  and  P'  for  which 
the  process  is  allowed.  The  total  observed  intensity  is 
then  equal  to  the  sum  of  all  such  y-ray  contributions. 

The  distorted  theoretical  angular  distribution  is  then 
given  by 


Wdis (9 ' Ji' Jf /6)  “ 


I  a  (E,0  )  exD(-y  Z  )  a  , (E  ) 
n  n  ^  Hn  n  nn '  n 


PP' 


(18) 


X  wth (En'^y ' Jif Jf /6)  exp(-y  £  )  dEdVdA 


where  an^E/®n^ 


n 


a  ,  (E  ) 
nn1  n' 


Y 


Wth(En'^y'Ji' Jf 


is  the  differential  cross  section  for 

the  neutron  production  reaction  at  an 

angle  0^,  for  incident  proton  energy  E. 

is  a  normalizing  function. 

is  the  neutron  attenuation  coefficient. 

is  the  total  cross  section  for  the  (n,n') 

reaction  for  neutrons  of  energy  E  . 

n 

is  the  y-ray  attenuation  coefficient. 

6)  is  the  theoretical  angular  distribution 
of  gamma  rays  from  the  -*■  Jf  transition, 
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having  a  mixing  ratio  6 ,  for  incident 
neutron  energy  E  ,  where  the  outgoing 
y-ray  makes  an  angle  of  with  respect 
to  the  incident  neutron  direction. 
dEdVdA  are  the  integration  variables  over  the 

proton  energy ,  the  scatterer  volume,  and 
the  area  of  the  beamspot. 

3  3 

The  differential  cross  sections  for  the  H(p,n)  He 
reaction  used  were  those  of  Perry  et  at.  (Pe60) .  The 
neutron  and  y-ray  attentuation  coefficients  were  taken 
from  standard  compilations  of  neutron  and  photon  cross 
sections  (Go56,  St67).  The  total  (n,n')  inelastic  cross 
section  was  taken  to  be  proportional  to  the  square  root 
of  the  difference  between  the  neutron  energy  and  the 
threshold  energy  (B152). 

The  theoretical  angular  distributions  were  calcul¬ 
ated  by  using  a  series  of  computer  programs.  First,  the 
transmission  coefficients  for  the  , 116cd  and  ^^Nd 
nuclei  were  calculated  by  the  computer  program  TRANCO, 
written  by  H.  S.  Sherif  (Sh76) .  The  transmission 
coefficients  are  given  in  terms  of  the  optical  model 
phase  shifts  60 .  by  the  relation 

T. .  =  1  -  | exp (2£5 n  .) | 2  (19 

Si  3  1  Slj 

The  phase  shifts  were  determined  by  matching  the  internal 
and  external  wave  functions  at  a  point  where  the  nuclear 


. 


26 


potential  is  zero.  The  potential  used  was  the  optical 

model,  with  the  parameters  of  Wilmore  and  Hodgson  (Wi64). 

The  transmission  coefficients  generated  by  TRANCO 

were  used  as  input  by  the  program  AK,  written  by 

P.  W.  Green  (Gr71) .  This  program  calculated  the  relative 

population  parameters  of  the  magnetic  substates  w (M  ) 

from  the  Compound  Nucleus  Statistical  Model.  Tables  of 

B  coefficients,  as  a  function  of  neutron  energy  and  the 
K 

initial  state's  spin  and  parity,  were  generated  for  each 

level  studied.  A  test  was  carried  out  to  determine  the 

dependence  of  B  coefficients  on  the  magnitude  of  the 

K 

transmission  coefficients.  In  this  test,  the  optical 

model  parameters  of  Wilmore  and  Hodgson  (Wi64),  Perey 

and  Buck  (Pe62) ,  Rosen  et  at.  (Ro65) ,  and  Bjorklund  and 

Fernbach  (Bj58)  were  used  to  calculate  transmission 

114 

coefficients  for  levels  m  Cd .  It  was  found  that, 

whereas  the  transmission  coefficients  generated  by  the 

different  sets  of  optical  model  parameters  differed  by 

as  much  as  70%,  the  angular  distribution  coefficients 

varied  by  only  about  3%.  These  results  are  in  agreement 

with  those  of  Pilte.;t  at  (Pi70)  ,  who  carried  out  a 

5  5 

similar  test  for  the  Mn(p,n)  reaction. 

The  B  coefficients  were  used  as  input  for  the 
K 

computer  program  EVA,  written  by  J.  M.  Davidson  (Da76b). 
EVA  numerically  integrated  Eq.  (18)  to  find  the  distorted 
theoretical  angular  distributions  for  specified  values 
J ,  and  6.  These  were  then  compared  to  the 


of  J  ,  tt  , 


experimental  data  by  means  of  the  x2  test.  Plots  of  x2 
as  a  function  of  the  multipole  mixing  ratio  (-90°  < 
arctan6  <  90°)  were  generated  for  all  possible  values  of 
the  initial  state's  spin-parity  consistent  with  dipole 
and  quadrupole  radiation.  The  0.1%  confidence  level 
(En74)  was  taken  as  the  criterion  for  rejecting  a  spin 
assignment . 

3 . 2  Experimental  Setup  and  Apparatus 

The  experimental  setup  is  shown  in  fig.  3.  Protons 
were  used  to  produce  neutrons  via  the  H(p,n)  He  reaction. 
The  tritium  targets  consisted  of  a  3.27  mg/cm2  layer  of 
metallic  erbium  deposited  on  a  tantalum  backing,  with 
tritium  absorbed  in  the  erbium  layer  at  an  atomic  ratio 
of  1:1.  A  10  nsec  pulsed  beam  with  a  repitition  rate  of 
1  MHz  and  an  average  current  of  .7  pA  was  provided  by 
the  University  of  Alberta  Van  de  Graaff  accelerator. 
Time-of-f light  gating  of  the  signals  from  each  of  two 
Ge(Li)  detectors  reduced  the  background  from  the  (n,n'y) 
reaction  in  the  detectors  and  gave  an  upper  limit  of 
approximately  20  nsec  for  the  lifetimes  of  the  observed 
states.  The  accumulated  charge  per  angle  for  the  six 
hour  runs  was  typically  15  mC. 

A  tantalum  collimator  with  an  aperture  of  0.24 
cm  diameter  was  placed  approximately  1  cm  from  the  neutron 
production  target.  The  cylindrical  scatterer  was  placed 
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Ge(Li) 


NE213 


Figure  3.  Experimental  set-up  and  electronics 
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a  distance  of  0.5  cm  from  the  neutron  production  target. 

It's  axis  and  the  beam  axis' were  colinear.  Angular 

distributions  were  recorded  for  all  nuclei  at  two 

different  proton  energies.  These  energies,  along  with 

the  scatterer  characteristics,  are  listed  in  Table  IV. 

Two  Ge(Li)  detectors  placed  about  50  cm  from  the 

scatterer,  detected  the  y-ray  radiation.  The  angles  0^ 

at  which  data  was  taken  were  90°,  65°,  50°,  40°,  25°, 

and  0°.  Neutron  flux  was  monitored  using  a  stationary 

NE213  liquid  scintillation  detector  placed  at  a  distance 

137 

of  several  meters  from  the  target.  A  Cs  source, 

=  661  keV,  was  placed  above  the  Ge(Li)  detectors  to 
provide  a  means  of  calibrating  dead  time  and  gain 
stability  of  the  system.  Signals  from  the  detectors 
were  analyzed  and  stored  on  line  using  a  Honeywell 
DDP-516  computer. 


Table  IV.  Scatterer  Characteristics. 
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CHAPTER  4 


RESULTS 


The  properties  of  the  low-lying  states  of  the 

114,116  ,  14  6  ,  . 

Cd  and  Nd  nuclei  were  determined  by  measuring 

the  angular  distributions  of  y-ray  transitions.  The 
experimental  data  was  then  compared  to  the  theoretical 
angular  distributions,  which  were  generated  for  all 
values  of  the  initial  state's  spin  and  parity  such  that 
the  y-ray  transition  was  consistent  with  dipole  and 
quadrupole  radiation.  The  sensitivity  of  the  angular 
distribution  to  parity  is  due  to  the  dependence  of  the 
coefficients  on  the  initial  state's  alignment.  For 
all  transitions  studied  in  this  work,  calculations  were 
done  for  both  parity  assignments.  In  no  case  was  it 
possible  to  exclude  one  parity  assignment  on  the  basis 
of  the  0.1%  confidence  level.  Since  the  differences 
are  usually  small,  all  results  shown  in  this  section, 
except  where  noted,  are  for  positive  parity  assignments 
Parity  assignments  were  made  on  the  basis  of  the 

+  1  per  degree  of  freedom  rule  and  a  consideration 


of  lifetimes 


4.1  Results  for 


114 

A  typical  y-ray  spectrum  for  Cd,  taken  at  a  pro¬ 
ton  bombarding  energy  of  3.60  Mev,  is  shown  in  fig.  4a. 
The  numbers  above  each  peak  denote  the  energies,  in  keV, 
of  the  y-ray  transitions  between  levels  in  114cd.  The 
deduced  level  and  decay  scheme,  as  shown  in  fig.  5,  is 
in  good  agreement  with  the  level  and  decay  scheme  pro¬ 
posed  by  Gill  et  al .  (Gi74).  Table  V  summarizes  the 

mixing  ratios  and  spin-parity  assignments  determined  in 

this  work.  The  errors  on  the  mixing  ratios  were  deter- 

2 

mined  by  the  x  •  +  1  rule  (Ro75) .  As  the  first  excited 

mm 

114 

state  of  cd  has  been  extensively  investigated  in 
previous  works,  it  has  not  been  studied  here. 

4.1.1  The  1134,  1210,  and  1283  keV  levels. 

These  three  levels  have  previously  been  shown  in 

several  experiments  (Ki75)  to  be  the  triplet  of  two 

phonon  states,  having  spin-parities  of  0  ,  2  ,  and  4  . 

The  results  for  the  1134  ->  558  keV  transition  are  shown 

2 

m  fig.  6.  As  the  plot  of  x  versus  arctanS  shows, 
the  results  are  consistent  with  a  spin-parity  assignment 
of  0+  for  the  1134  keV  state.  The  level  at  1210  keV 
was  seen  to  decay  via  y-rays  of  651  and  1210  keV  to  the 
first  excited  and  ground  states.  The  1210  ->  0  keV 
transition,  however,  was  seen  in  sufficient  strength 
only  during  the  Er  <  2.67  MeV  runs,  when  this  state  was 
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Figure  4. 


Typical  y-ray  spectra  for  the  Cd  isotopes,  a)  114Cd.  The  numbers 
above  each  peak  refer  to  the  energies,  in  keV,  of  the  y-ray  tran¬ 
sitions  between  the  Cd  levels.  Background  peaks  are  denoted  by  B. 
b)  116Cd.  The  inset  shows  the  y-ray  decays  of  the  triplet  of  states 
near  1920  keV. 
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Figure  5.  The  deduced  level  and  decay  scheme  for  4Cd.  The 

dotted  lines  refer  to  transitions  that  were  too  weak 
to  provide  useful  information.  The  spin-parity  assign¬ 
ments  are  the  result  of  the  present  and  previous 
works . 
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Table  V.  The  Gamma 

Ray  Mixing 

Ratios  and 

Spin-Parity 

Assignments  for  "^^Cd 

Determined 

in  This  Work. 

,  a) 

Level  Energy  y- 

ray  Energy 

j1!  +  jI 

i  f 

ArctanS 

1134 

576 

0+  -  2+ 

0 

1210 

651 

2+  ■*  2  + 

58  13 

1283 

725 

4+  ->  2+ 

ob) 

1365 

806 

2+  2+ 

-80  5 

16  5 

1733 

1175 

4+  -*  2+ 

0b) 

552 

4+  -  2+ 

ob> 

1843 

1285 

2+  -  2+ 

°l 

1 — 1 1 

O 

00 

1 

21  15 

708 

2  ->  0 

0 

1861 

1303 

3+  2+ 

62  ->  16 

1865 

1307 

0  -y  2+ 

0 

1  +  2  + 

all  values 

2  -y  2  + 

5  ->  90 

3  -+  2  + 

-15  ->  20 

55  +  90 

4+  +  2+ 

0b) 

1959 

1401 

1+  -►  2+ 

-50  35 

2+  +  2+ 

52  14. 

3+  -v  2+ 

72  12 

(continued . . . . ) 
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Table  V  (continued) 
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a)  The  energies  are  in  keV.  The  uncertainties  are  ±1  keV. 

b)  The  radiation  was  assumed  to  be  pure  quadrupole. 


Arctan  6 

Figure  6.  The  1134  keV  level:  a)  angular  distribution  for  the 
576  keV  y-ray,  and  b)  x2  curves. 
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fed  through  y-ray  cascades  from  higher  excited  states. 

The  results  for  the  1210  -*  558  keV  transition,  shown  in 
fig.  7,  are  consistent  with  a  spin-parity  assignment  of 
2+  for  this  level.  The  angular  distribution  measurements 
for  the  1283  -*  558  keV  transition,  shown  in  fig.  8,  are 
consistent  with  an  initial  state  spin  parity  of  4+. 

4.1.2  The  1306  and  1365  keV  levels. 

The  nature  of  these  states  is  not  very  well  known, 
but  it  has  been  proposed  that  they  may  be  considered  as 
quasi  (single)  particle  states  (Sa64a) .  The  1306  keV 
level  has  two  reported  y-ray  decay  modes,  to  the  1210 
and  558  keV  levels  (Ba66)  .  The  1306  -*■  1210  keV  transi¬ 
tion  was  not  observed  as  the  y-ray  energy  spectrum  had  a 
cutoff  of  approximately  250  keV.  The  1306  ->  558  keV 
transition  was  observed  in  this  work,  but  was  too  weak 
to  provide  useful  information.  The  weakness  of  this 
transition  is  probably  due  to  competition  from  the  decay 
to  the  ground  state  via  internal  conversion,  as  observed 
by  Backlin  et  al .  (Ba66) .  The  1365  ■*  558  keV  transition 

was  observed  in  this  work.  The  angular  distribution 
measurements,  shown  in  fig.  9,  are  consistent  with  a 
spin-parity  assignment  of  2+  that  has  been  suggested  by 
earlier  work  (Ki75) . 


GAMMA  RAY  COUNTS 


b) 


Figure  7 


T 


The  1210  keV  level:  a)  angular  distribution  for  the 
651  keV  y-ray,  and  b)  x2  curves. 


a) 


Figure  8.  The  1283  keV  level:  a)  angular  distribution  for 
725  keV  y-ray,  and  b)  x2  curves. 


a) 


5500 


FigurG  9..  ThG  1365  keV  level i  a)  ancjular  distribution  for  thfi 
806  keV  Y-ray,  and  b)  x2  curves. 
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4.1.3  The  1733,  1843/  1861,  and  1865  keV  levels. 

These  states  all  lie  well  above  the  triplet  of 
two-phonon  quadrupole  states,  and  well  below  the  quad- 
rupole-octupole  states  suggested  by  Sakai  and  Tamura 
(Sa64a) .  Curiously,  these  states  were  not  observed  in 
in  Coulomb  excitation,  which  should  preferentially  ex¬ 
cite  collective  states. 

The  state  at  1733  keV  was  seen  to  decay  to  the 
levels  at  558,  1210,  and  1306  keV  excitation.  The  angular 
distribution  measurements  for  the  two  stronger  transitions, 
shown  in  figs.  10  and  11,  placed  no  restrictions  on  the 
value  of  the  spin  of  this  state.  A  spin  assignment  of 
4+  has  previously  been  made  (Gr68,Gi74). 

Four  decay  modes  have  been  reported  for  the  level 
at  1843  keV  excitation  (Ba66).  Of  these,  only  the  transi¬ 
tions  to  the  levels  at  1134  and  558  keV  were  seen.  As 
the  1134  keV  state  has  a  spin  of  0+,  the  multipole  mixing 
ratio  for  the  1843  ■+  1134  keV  transition  must  be  zero. 

As  fig.  12  shows,  the  only  spin  for  which  the  value  of 
X  (6  =  0)  was  below  the  0.1%  confidence  limit  was  2.  If 
this  state  has  negative  parity,  the  transition  to  the 
558  keV  level,  a  2+  state,  would  be  almost  pure  El, 

giving  a  6  of  zero.  The  deduced  values  of  arctan6  of 

2 

±(80°  ±  10°)  and  21°  ±  12°  from  the  plot  of  x  versus 
arctanS,  shown  in  fig.  13,  make  this  assignment  relatively 
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Figure  10.  The  1733  keV  level:  a)  angular  distribution  for  the 
1175  keV  y-ray,  and  b)  x*  curves. 
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a) 
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b) 


Figure  11.  The  1733  keV  level:  a)  angular  distribution  for  the 
522  keV  y-ray,  and  b)  x2  curves. 


GAMMA  RAY  COUNTS 


Figure  12.  The  1843  keV  level:  a)  angular  distribution  for  the 
708  keV  y-ray,  and  b)  x2(<$=0). 
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Figure  13.  The  1843  keV  level:  a)  angular  distribution  for  the 
1285  keV  y-ray,  and  b)  x2  curves. 
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improbable.  Therefore,  we  assign  a  spin-parity  of  2+ 
to  the  1843  keV  state. 

The  results  for  the  1861  -*■  558  keV  transition  are 

2 

shown  in  fig.  14.  As  this  figure  shows,  minima  in  x  (5) 
permit  only  an  initial  state  spin  of  3.  The  values  of 
arctan  6  of  16°  -*■  62°  obtained  from  the  )^n  +  1  rule 

make  it  relatively  improbable  that  this  is  the  octupole 
state.  The  spin  assignment  of  3+  suggests  that  the 
1861  keV  state  is  one  of  the  three-phonon  quadrupole 
states.  This  state  has  been  observed  in  the  study  of 
Cd(n,n'y)  excitation  functions  (Gi74) ,  and  was  tentatively 
assigned  a  spin  of  1. 

The  excited  state  at  1865  keV  has  one  reported 
decay  mode;  to  the  first  excited  state.  The  angular 
distribution,  shown  in  fig.  15,  was  found  to  be  isotropic, 
which  placed  no  restriction  on  the  spin  of  this  state. 

4.1.4  The  1959,  2049,  and  2206  keV  levels. 

The  1959  keV  state  has  reported  y-ray  branches  to  the 

established  2+  states  at  558,  1210,  and  1365  keV.  The  594 

keV  y-ray  from  the  1959  -*■  1365  keV  transition  could  not  be 

74 

resolved  due  to  a  background  peak  from  the  Ge(n,n  y) 
reaction.  This  was  also  the  case  in  the  (n,n'y)  excitation 
work  of  Gill  zt  a£.(4)  Although  the  748  keV  y-ray  has  also 
been  placed  in  the  decay  scheme  as  the  1306  -*  558  keV  tran¬ 
sition,  it  has  previously  been  shown  (4)  that  the  1959  1210 


49 


Figure  14.  The  1861  keV  level:  a)  angular  distribution  for  the 
1303  keV  y-ray,  and  b)  x2  curves. 


50 
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Figure  15.  The  1865  keV  level:  a)  angular  distribution  for  the 
1307  keV  y-ray,  and  b)  x2  curves. 


keV  transition  is  mainly  responsible  for  the  observed 

strength.  Sample  angular  distributions  for  the  1959  -*■ 

1210  and  1959  558  keV  transitions  and  plots  of 

versus  arctan  6  are  shown  in  Fig.  6.  As  this  Figure 

2 

shows,  minima  in  x  occur  below  the  0.1%  confidence 

level  only  for  initial  spins  of  1,2,3  for  both  transi- 

2 

tions.  The  minimum  x  per  degree  of  freedom  for  the 
spin  assignment  1  corresponds  to  a  confidence  level  of 
only  0.5%,  making  that  assignment  relatively  improbable. 

It  should  also  be  noted  that  for  a  spin  of  3,  the  deduced 
values  of  arctan  6  are  consistent  with  a  negative  parity 
assignment . 

The  1959  keV  state  was  originally  assigned  a  spin- 
parity  of  3  in  the  Coulomb  excitation  work  of  Sakai  Q,t 
al.(7)  This  assignment  was  later  disputed  by  the  (n,n'y) 
work  of  Gill  at  aZ  (4)  whose  measured  cross  section  for 
the  1959  keV  state  was  too  low  to  support  a  3"  assignment. 
The  fact  that  the  594  keV  y-ray  could  not  be  resolved  made 
the  assignment  of  J  =  1  in  that  work  for  the  1959  keV 
level  tentative  only.  This,  taken  in  conjunction  with  the 
0.5%  confidence  level  for  a  spin  of  1,  and  the  fact  that 
the  1959  keV  state  was  observed  through  Coulomb  excitation 
a  reaction  that  should  favor  a  single  phonon  transition, 
leads  us  to  assign  a  spin-parity  of  3  . 

The  y-ray  decay  of  the  2049  keV  state  to  the  first  ex 
cited  state  was  observed  in  this  work.  The  angular  distri 
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b) 


1959  *  1210  keV 


Figure  16.  The  1959  keV  level:  a)  angular  distribution  for  the 
748  keV  y-ray,  and  b)  x2  curves. 
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Figure  17.  The  1959  keV  level:  a)  angular  distribution  for 
1401  keV  y-ray,  and  b)  x2  curves. 
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Figure  18.  The  2049  keV  level:  a)  angular  distribution  for  the 
1491  keV  y-ray,  and  b)  x2  curves. 
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bution  measurements,  shown  in  fig.  18,  allowed  an  initial 
state  spin  of  3  with  an  arctan  of  the  multipole  mixing 
ratio  6  of  70°  ±  10°,  implying  a  spin-parity  of  3+  for 
this  level.  Thus  the  assignment  is  3+.  The  1648  keV 
Y*-ray  from  the  2206  ->  558  keV  transition  was  too  weak  to 
provide  any  useful  information. 

4 . 2  Results  for  llecd. 

116 

A  typical  y-ray  spectrum  for  Cd  is  shown  in  fig. 

4b.  The  deduced  level  and  decay  scheme  is  shown  in  fig. 

114 

19.  As  in  the  case  of  Cd,  the  first  excited  state  of 
116 

Cd  was  not  studied  in  this  work.  The  values  of  the 
multipole  mixing  ratios  and  the  spin-parity  assignments 
deduced  in  this  work  are  summarized  in  Table  VI. 

4.2.1.  The  1213,  1220,  1283,  and  1381  keV  levels. 

The  first  two  of  the  two-phonon  quadrupole  states  are 
at  1213  and  1220  keV  excitation  respectively.  The  primary 
decay  mode  for  these  states  is  by  y-ray  de-excitation  to 
the  level  at  513  keV.  The  observed  branching  ratios  for 
the  1213  keV  level  are  56  ±  7%  for  the  1213  513  keV 

transition  and  44  ±  7%  for  the  1213  -*■  0  keV  transition. 

The  error  for  these  results  is  associated  with  the  non¬ 
isotropic  weighting  of  the  angular  distribution.  The  1220 
keV  state  was  observed  to  decay  only  to  the  513  keV  state. 
The  angular  distribution  measurements  for  the  1213 
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19.  The  deduced  level  and  decay  scheme  of  116Cd.  The 

dotted  lines  refer  to  those  transitions  that  were  ob¬ 
served,  but  were  too  weak  to  provide  useful  infor¬ 
mation.  The  spin-parity  assignments  are  the  result 
of  the  present  and  previous  works. 
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Table  VI.  The  Gamma  Ray  Mixing  Ratios  and  Spin-Parity 

1  1 

Assignments  for  cd  Determined  in  This  Work. 


^  \ 

Level  Energy 

y-ray  Energy 

jI  +  jj 

Arctan6 

1213 

1213 

2+  ->  0+ 

0 

700 

2+  ^  2+ 

57  13 

1220 

706 

4+  -  2+ 

ob) 

1283 

769 

0+  ->  2  + 

0 

l±  -  2  + 

all  values 

2+  -  2+ 

28  ->  90 

3+  -►  2+ 

78  8 

3"  +  2  + 

-2  7 

1381 

867 

0+  ->  2+ 

0 

1644 

1644 

2+  -  0+ 

0 

1130 

2+  -  2+ 

-32  19 

1917 

1404 

3+  -  2+ 

38  2_0 

1923 

1409 

1+  ->  2+ 

-51  2_5 

2+  ->  2  + 

63  18 

3+  -  2  + 

75  9 

3~  ■>  2+ 

8  8 

1930 

1416 

0+  -  2+ 

0 

1+  ->  2+  all  values 

2±  ->  2+  -90  ->  -52;  -7  -*  90 

3±  ->  2+  -11  16 


(continued . . . . ) 
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Table  VI  (continued) 


2120 


1607 


4+  -  2+ 
0+  *  2+ 


ob) 

0 

4  +  69 


2  + 


3" 


4 


+ 


1  ii 

-19  8 


a) ,  b)  The  notation  is  as  in  Table  V. 
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Figure  20.  The  1213  keV  level:  a)  angular  distribution  for  the 
700  keV  y-ray,  and  b)  x2  curves. 
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Figure  21.  The  1213  keV  level:  a)  angular  distribution  for  the 
1213  keV  y-ray,  and  b)  x2(<$). 
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513  and  1213  -*  0  keV  transitions  are  shown  in  figs.  20 
and  21  respectively.  These  results  are  consistent  with 
an  initial  state  spin  and  parity  of  2+  for  the  level  at 
1213  keV  excitation.  The  data  for  the  1220  -*■  513  keV 
transition  is  shown  in  fig.  22.  The  previous  assignment 
(Ca75)  of  4+  to  the  1220  keV  state  is  consistent  with  the 
results  of  the  present  work. 

Previous  information  about  the  level  at  1283  keV  ex¬ 
citation  comes  from  the  (p,p')  | ref .  De72 | ,  (p,p'y)  |ref. 

Be73|,  and  (n,n'y)  | ref .  Gi74 |  reactions.  Spin  assign¬ 

ments  made  on  the  basis  of  these  experiments  were  2+ ,  2+, 
and  0+  respectively.  As  fig.  23  shows,  the  angular  dis¬ 
tribution  for  the  1283  -*  513  keV  transition  is  consistent 
with  initial  state  spins  of  0,  1,  2,  and  3. 

Angular  distribution  measurements  for  the  1381  -*■  513 
keV  transition  are  consistent  with  a  spin  assignment  of 
0+  for  the  level  at  1381  keV  excitation.  The  angular  dis- 

7T  + 

tribution  is  shown  in  fig.  24.  An  assignment  of  J  =0 
for  this  level  has  been  made  in  previous  works  (Ca73) . 

4.2.2  The  1644  keV  level. 

In  the  present  work,  we  have  observed  a  level  at  1644 
keV  de-excited  by  y-rays  of  1644  and  1130  keV  to  the 
ground  and  first  excited  states  respectively.  The  results 
for  these  transitions  are  shown  in  figs.  25  and  26.  The 
1644  -v  0  keV  transition  allows  possible  spin  assignments 
for  this  state  of  2  and  3.  Values  of  arctan  6  for  the 


‘ 


GAMMA  RAY  COUNTS 


62 


Figure  22.  The  1220  keV  level:  a)  angular  distribution  for  the 
706  keV  y-ray,  and  b)  x2  curves. 
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Figure  23.  The  1283  keV  level:  a)  angular  distribution  for  the 
769  keV  y-ray,  and  b)  X2  curves. 
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Figure  24.  The  1381  keV  level:  a)  angular  distribution  for  the 
867  keV  y-ray,  and  b)  x2  curves. 
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Figure  25.  The  1644  keV  level:  a)  angular  distribution  for  the 
1644  keV  y-ray,  and  b)  y2  curves. 
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Figure  26.  The  1644  keV  level:  a)  angular  distribution  for  the 
1130  keV  y-ray,  and  b)  y2  curves. 
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1644  -►  513  keV  transition  (J^  -*■  2+)  further  limit 
the  initial  state  spin-parity  to  2+  or  3  +  .  If  this 
state  is  3+,  then  the  decay  of  the  1644  keV  state  to 
the  ground  state  will  be  by  M3  radiation  only.  This  is 
relatively  improbable  as  the  Weisskopf  estimate  for  the 
lifetime  of  this  state  is  5  ysec.  Therefore,  we  deduce 
that  the  1644  keV  level  has  a  spin  of  2+.  A  tentative 
assignment  of  (2,3)  has  previously  been  made  (Gi74). 

4.2.3  The  1917,  1923,  and  1930  keV  levels. 

The  states  at  1917,  1923,  and  1930  keV  excitation 

occur,  within  errors,  where  the  3_  octupole  state  was 

originally  identified  [ref.  Ko69,  De72].  These  states 

all  have  reported  (Gi74)  y-ray  branches  to  the  513  keV 

first  excited  state,  which  has  an  established  spin-parity 

of  2+.  Sample  angular  distribution  data  and  plots  of 
2 

X  (6)  for  these  three  y-ray  transitions  are  shown  in 

Figs.  27,  28,  and  29,  respectively.  For  the  1917  513 

keV  transition,  minima  occur  below  the  0.1%  confidence 

level  only  for  an  initial  state  spin  of  3,  in  disagreement 

with  a  previous  (Gi74)  tentative  assignment  of  0.  The 

large  quadrupole  strength  (arctan6  =  38°  ±  20°)  for  this 

transition  makes  it  unlikely  that  the  1917  keV  state  is 

2 

3  .  For  the  1923  -*  513  keV  transition,  minima  in  x  (6) 
occur  below  the  0.1%  confidence  limit  for  initial  spins 
of  1,  2,  and  3.  The  value  of  arctan6  of  -8°  ->  8°  deduced 
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Arctan  6 

Figure  27.  The  1917  keV  level:  a)angular  distribution  for  the 
1404  keV  y-ray,  and  b)  the  x2  curves. 
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Figure  28.  The  1923  keV  level:  a)  angular  distribution  for  the 
1409  keV  y-ray,  and  b)  x2  curves. 


1930  — 


Figure  29.  The  1930  keV  level:  a)  angular  distribution  for 
1416  key  y-ray,  and  b)  x2  curves. 
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for  an  initial  spin  of  3  is  consistent  with  an  assignment 
of  3  for  the  1923  keV  state.  The  data  for  the  1930  -> 

513  keV  transition  place  no  restriction  on  its  spin- 
parity  . 

4.2.4  The  1953  ,  2044  ,  and  2120  keV  levels. 

The  1953  -*  513  keV  transition  was  observed,  but 
the  large  background  in  the  forward  angles  prohibited 
any  useful  information  from  being  deduced.  The  2044 
513  keV  transition  was  too  weak  to  provide  any  useful 
information.  The  angular  distribution  measurements, 
for  the  2120  -»•  513  keV  transition,  which  are  shown  in 
fig.  30,  placed  no  restriction  on  the  spin-parity  of 
the  2120  keV  state. 

146 

4 . 3  Results  for _ Nd . 

146 

Typical  y-ray  spectra  for  the  Nd  experiment  are 
shown  in  fig.  31.  The  y-ray  transitions  between  the 
levels  in  Nd  are  labeled  by  their  energies,  in  keV. 

The  energies  determined  in  this  work,  with  the  exception 
of  the  level  at  1377  keV,  previously  given  as  1372  keV, 
are  in  excellent  agreement  with  those  adopted  by 
Burrows  (Bu75) .  The  level  and  decay  scheme  is  shown  in 
fig.  32.  This  scheme  differs  from  that  of  Burrows  in 
that  a  new  level  is  proposed  at  1787  keV,  while  the 
previously  adopted  levels  at  1515,  1538,  1751,  1837, 
and  1913  keV  are  absent  as  they  were  not  observed  in  this 
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Figure  30.  The  2120  keV  level:  a)  angular  distribution  for  the 
1607  keV  y-ray,  and  b)  x2  curves. 
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Figure  31.  Typical  y-ray  spectra  for  the  146Nd  experiment, 
taken  at  En^2.67  MeV,  0=90°.  a)  The  scatterer 
in  spectrum.  The  numbers  above  each  peak  refer 
to  the  energies,  in  keV,  of  the  y-ray  transitions 
in  146Nd.  The  three  peaks  that  have  lines  on 
each  side  are  drawn  at  1/2  of  the  observed  strength, 
b)  The  scatterer  out  spectrum.  The  data  collec¬ 
tion  time  for  this  spectrum  was  approximately 
1/2  of  the  time  taken  for  the  spectrum  shown  in  a). 


Figure  32. 


The  deduced  level  and  decay  scheme  of  146Nd.  The 
dotted  lines  refer  to  transitions  that  were  ob¬ 
served,  but  were  too  weak  to  provide  useful  infor¬ 
mation.  The  spin-parity  assignments  are  the  re- 
sv 1 t  of  the  present  and  previous  works. 
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work.  The  spin-parity  assignments  and  multipole  mixing 
ratios  deduced  in  this  work  are  summarized  in  Table  VII. 

The  one  phonon  2+  state  at  454  keV  was  not  studied  in 
this  work. 

4.3.1  The  998 ,  1043 ,  and  1049  keV  levels. 

These  three  states  lie  at  approximately  twice  the 
energy  of  the  one-phonon  2+  state.  Therefore,  one 
might  expect  that  they  are  members  of  the  two-phonon 
triplet.  The  998  keV  has  two  reported  decay  modes;  to 
the  ground  and  first  excited  states  (Gr68) .  Although 
neither  of  these  transitions  was  seen  in  this  work,  a 
gamma  ray  transition  from  the  state  at  1471  keV  excita¬ 
tion  to  this  level  was  observed.  The  angular  distribution 
for  this  transition  is  shown  in  fig.  33.  This  data  was 
put  through  the  computer  program  EVA  with  fixed  at 

2  and  allowed  to  vary  from  0  to  4 .  The  resultant 
2 

X  (6)  plot  allowed  possible  spins  of  1,  2,  and  3  for 
this  state.  The  results  for  the  1043  454  keV  trans¬ 

ition,  shown  in  fig.  34,  are  consistent  with  the  previous 
assignment  (Bu75)  of  4+  for  the  level  at  1043  keV 
excitation.  The  1049  keV  level  has  one  reported  decay 
mode;  to  the  454  keV  level  (Gr68).  Unfortunately,  the 

large  background  in  the  y-ray  spectrum  due  to  the 
74 

Ge(n,n'y)  reaction  prohibited  any  useful  information 


from  being  deduced. 
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Table  VII.  The  Gamma  Ray  Mixing  Ratios  and  Spin-Parity 

146 

Assignments  for  Nd  Determined  in  This  Work. 


Level  Energy 


a) 


Y-ray  energy 


1 


-> 


.TT 


ArctanS 


998 


1043 

1189 

1377 

1471 


474 


590 

735 

1377 

924 

1471 


2+  -  1 
2+  +  2 
2+  +  3 
4+  ->  2+ 


3  ->2 


-»■ 


0 


+ 


1  ->  2 


+ 


0 


+ 


10  +  60 
37  -  73 


-11  +  -23 


0 


b) 


5  I 
0 


-90  +  90 
0 


1777 


1778 


1978 


1323 


1333 


1524 


0+  0+ 

2  2 


-32  18 


3+  ■>  2+ 
4+  ->  2+ 
0+  ->  2+ 


1  -> 

2  -> 


3  2 


0 


+ 


+ 


4+  2+ 


,+ 


1+  ->  2+ 


2  ->  2 


+ 


-81  +  -25 

ob) 

0 

-90  ->  90 
-81  ->  27 
-16  11 
ob) 

0 

±  (90  ■>  45) 
-75  -*  14 


(continued . . . . ) 
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Table  VII  (continued) . 
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a) ,  b)  notation  is  as  in  Table  V. 
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Figure  33.  The  999  keV  level:  a)  angular  distribution  for  the 
474  keV  y-ray,  and  b)  y2  curves. 
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Figure  34.  The  1043  keV  level:  a)  angular  distribution  for  the 
590  keV  y-ray,  and  b)  x2  curves. 
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4 • 3 • 2 ♦  The  1189 ,  1377,  1471,  and  1687  keV  levels. 

The  data  and  results  for  the  1189  -*  454  keV  transi¬ 
tion  are  shown  in  fig.  35.  As  the  x^(6)  plot  illustrates, 
the  results  are  consistent  with  the  previous  spin-parity 
assignment  of  3  (Bu75) .  Therefore,  for  ^^Nd,  the  oc- 

tupole  state  occurs  at  approximately  the  energy  predicted 
by  the  harmonic  vibrational  model.  The  level  at  1377  keV 
was  seen  to  decay  via  y-rays  of  1377  and  924  keV  to  the 
ground  and  first  excited  states  respectively.  The  angu¬ 
lar  distribution  measurements  for  the  1377  -*■  0  keV  transi¬ 
tion,  shown  in  fig.  36,  allowed  a  unique  spin  assignment 
of  1.  The  deduced  values  of  the  multipole  mixing  ratio 
for  the  1377  ->  454  keV  transition  (see  fig.  37)  did  not 
limit  the  parity  of  this  state.  The  level  at  1471  keV  ex¬ 
citation  was  seen  to  decay  via  y-rays  of  474,  1017,  and 
1471  keV.  The  474  keV  y-ray  has  already  been  discussed  in 
the  context  of  the  999  keV  level.  In  the  forward  angle 

runs,  we  were  not  able  to  resolve  the  1017  keV  y-ray  from 
27 

the  A1  background  peak.  The  results  for  the  1471  0 

keV  transition,  shown  in  fig.  38,  permit  initial  state 
spins  of  2  or  3.  A  spin  assignment  of  2+  has  previously 
been  suggested  (Bu75) .  The  level  at  1687  keV  was  observed 
to  decay  to  the  one  phonon  2+  state  by  a  y-ray  of  1233  keV. 
As  can  be  seen  in  the  scatterer  out  y-ray  spectrum  (see 
fig.  31b)  this  is  the  same  energy  as  a  y-ray  transition  in 
19F.  The  statistics  for  the  146Nd  1233  keV  y-ray  that 
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Figure  35.  The  1189  keV  level:  a)  angular  distribution  for  the 
735  keV  y-ray,  and  b)  x2  curves. 
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Figure  36.  The  1377  keV  level:  a)  angular  distribution  for  the 
1  377  keV  y-ray,  and  b)  x2(<$=0). 
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Figure  37.  The  1377  keV  level:  a)  angular  distribution  for  the 
924  keV  y-ray,  and  b)  x2  curves. 
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Figure  38.  The  1472  keV  level:  a)  angular  distribution  for  the 
1471  keV  and  b)  x2(d=0). 
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results  from  a  scatter-in  minus  scatterer-out  spec¬ 
trum  were  too  poor  to  allow  any  information  to  be 
deduced.  In  (n,y)  capture  work  (Gr68) ,  an  observed 
y-ray  of  1687  keV  was  attributed  to  this  level.  This 
was  later  contradicted  by  the  results  from  the  (d,d') 
reaction  (Ch70)  which  gave  a  spin-parity  assignment  of 
0+  to  the  1687  keV  level.  No  y-rays  of  energy  near 
1687  keV  were  observed  in  this  work. 

4.3.3  The  1777,  1787,  and  1978  keV  levels. 

The  results  for  the  1777  454  keV  transition  are 

shown  in  fig.  39.  The  values  of  the  initial  state's 
spin  that  are  below  the  0.1%  confidence  level  are  2,  2, 
and  4.  The  deduced  values  of  arctan6  for  these  spins 
further  limit  the  parity  to  be  positive.  No  spin  assign¬ 
ment  for  this  level  has  previously  been  made.  The  1333 
keV  y-ray,  as  with  the  1323  keV  y-ray  from  the  1777  -> 

454  keV  transition,  was  seen  only  during  the  E  =3.60 

P 

MeV  runs.  As  a  y-ray  of  this  energy  was  not  observed  in 

the  scatterer-out  spectrum,  the  1333  keV  y-ray  must  be 

the  result  of  a  transition  in  the  Nd  target.  No  y-rays 

of  this  energy  have  been  reported  for  the  other  Nd 

146 

isotopes.  Therefore,  a  new  level  for  Nd  has  been 
proposed  to  account  for  the  1333  keV  y-ray.  The  angular 
distribution  measurements  for  this  transition,  shown  in 
fig.  40,  place  no  restriction  on  the  spin  of  this  state. 
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Figure  39.  The  1777  keV  level:  a)  angular  distribution  for  the 
1323  keV  y-ray,  and  b)  x2  curves. 
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Figure  40.  The  1787  keV  level:  a)  angular  distribution  for  the 
1333  keV  y-ray,  and  b)  y2  curves. 


The  data  and  results  for  the  1978  454  keV  transi¬ 


tion  are  shown  in  fig.  41.  No  information  concerning 
the  spin  of  the  1978  keV  level  has  been  deduced. 
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.  The  1978  keV  level:  a)  angular  distribution  for  the 
1524  keV  y-ray,  and  b)  x2  curves. 


Figure  41 


CHAPTER  5 


DISCUSSION 

As  noted  in  the  introduction,  Scharff,  Goldhaber  and 
Weneser  (Sc55)  pointed  out  the  vibrational  nature  of  the 
first  and  second  excited  states  of  many  even-even  nuclei. 
It  is  my  intention  here  to  discuss  the  nature  of  the 
higher  excited  states  of  ^^'^^Cd  and  "^^Nd  with  respect 
to  the  predictions  of  the  harmonic  vibrational  model.* 
Considering  only  the  quadrupole  and  octupole  modes 
of  vibration,  the  spectrum  of  an  harmonic  vibrational 
nucleus  is  composed  of  highly  degenerate  and  equally 
spaced  levels.  The  levels  can  be  denoted  by 

|  X,N,JM  >  (2 

where  A  corresponds  to  the  mode  of  vibration,  N  the 

number  of  phonons,  J  the  total  angular  momentum,  and  M 

it's  z-pro j ection .  States  that  are  the  result  of  the 

coupling  of  phonons  of  different  modes  of  vibration  are 

denoted  by  |  A  =  A  +  X  ,N,JM  >.  The  levels  of  H4,116 
1  n  m 

146 

and  Nd  that  have  been  tentatively  characterized  as 


The  harmonic  vibrational  model  assumes  a  Hamiltonian  of 
the  form  H  =  T(a,  ;ir,  )  +  V(a.)  where  T  and  V  are  the 

A  A  A 

kinetic  and  potential  energies.  Only  second  order  terms 
of  a-y  and  are  used.  For  a  more  detailed  description 
of  the  vibrational  model,  see  Appendix  A. 
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vibrational  states  are  summarized  in  Table  VIII.  The 

bases  for  these  assignments  are  the  observed  energies 

and  deduced  spins  and  parities  of  the  levels  in  question. 

It  should  be  noted  that  these  three  nuclei  all 

possess  a  0+  ground  state  and  a  2+  first  excited  state. 

Also,  the  two  phonon  triplet  occurs  at  an  energy  close  to 

twice  the  energy  of  the  one  phonon  2+  state.  Taking  an 

average  of  the  triplet  energy  to  that  of  the  first 

excited  state  yields  the  ratios  of  2.17  for  ^^cd,  2.41 
116  146 

for  Cd,  and  2.29  for  Nd.  The  triplet  is  not 

14  6 

degenerate  and  in  Nd  two  of  the  members  of  the  triplet 
are  not  positively  identified,  while  in  the  Cd  isotopes, 
additional  states  occur  in  the  neighborhood  of  the 
triplet.  The  fact  that  members  of  the  triplet  are 
missing  may  be  due  to  experimental  shortcomings.  Two 
theories  have  been  put  forward  to  account  for  the  addi¬ 
tional  states.  Sakai  and  Tamura  (Sa64a)  postulated  that 
these  additional  states  may  be  thought  of  as  single 
particle  states  while  Eisenberg  and  Greiner  (Ei70) 
suggested  that  they  were  members  of  a  rotational  band  built 
on  a  deformed  excited  state. 

The  identification  of  the  three  phonon  states  is 
very  sparse  for  Nd.  For  Cd,  the  situation  is 

clearer,  but  the  |  2 , 3 , 6M  >  state  is  still  not  identified. 

TT  -f 

In  this  work,  a  state  with  J  =6  would  probably  not 
be  observed  as  the  transmission  coefficients  are  too 


93 


Table  VIII.  Vibration  Levels  of  , 116Cd  an£ 


| X ,  N,  JM  > 

114 

Cd  (keV) 

116Cd  (keV) 

146 

Nd  (keV) 

A 

O 

O 

o 

o 

0 

0 

0 

12,1,2m  > 

558 

513 

454 

|  2 , 2 , OM  > 

1134 

1283 

- 

| 2 , 2 , 2M  > 

1210 

1213 

- 

| 2,2,4M  > 

1283 

1220 

1043 

| 2,3,00  > 

(1865) 

(1930) 

- 

|  2 , 3 , 2M  > 

1843 

1644 

(1471) 

| 2 , 3 , 3M  > 

1861 

1917 

- 

| 2,3,4M  > 

1733 

- 

- 

|  2 , 3 , 6M  > 

- 

- 

- 

| 3,1,3M  > 

(1959) 

(1923) 

1189 

| 2+3, 2, JM  > 

— 

- 

(1377) 

Other  States 

1306  (0+) 

1381 

999 

Observed  in 

1365  (2+) 

1953 

1049 

This  Work 

2049  (3+) 

2044 

1777 

2246 

2120 

1787 

1978 

small  to  sufficiently  populate  the  level.  The  assignment 
of  the  1865  keV  and  1930  keV  levels  as  the  |  2 , 3 , 0M  > 
states  in  and  ^^Cd  respectively,  is  tentative  only 

For  cd  and  "^^Cd,  the  levels  at  1959  keV  and 
1923  keV  respectively  were  assigned  as  the  octupole  states 
As  the  multiplet  of  |  2+3, 2, JM  >  states  is  therefore 

expected  at  approximately  2500  keV,  they  were  not  observed 
This  was  not  the  case  for  °Nd  as  the  octupole  state 
was  seen  at  1189  keV.  The  state  at  1377  keV,  with  a  spin 
of  1,  is  therefore  at  approximately  the  energy  (454  +  1189 
keV)  that  one  would  expect  a  quadrupole-octupole  state. 

Nevertheless,  the  splitting  of  the  two  phonon  states 
and  also  their  deviation  from  the  predicted  energy 
spacing  are  indications  of  anharmonic  terms  (higher  order 
terms  in  a^)  in  the  Hamiltonian.  The  harmonic  approxima¬ 
tion  seems  to  be  a  zero-order  approximation  for  these 
nuclei . 

In  summary,  the  measurement  of  the  angular  distribu¬ 
tions  of  y-rays  from  the  H-4,116^^  ^4^Nd  (n,n'y)  reac¬ 
tions  enable  spin-parity  assignments  to  be  made  for  these 
nuclei  as  tabulated  in  Tables  V,  VI,  and  VII. 

In  114Cd,  the  states  at  1843  and  1861  keV  have  been 
assigned  spin-parities  of  2+  and  3+  respectively,  in  dis¬ 
agreement  with  previous  work  (6:74).  The  state  at  2049 
keV  has  been  assigned  a  spin  of  3+ .  While  the  state  at 
1959  keV  excitation  cannot  be  rigorously  limited  to  3  , 
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it  has  been  shown  that  the  previous  assignment  of  1  (6.74) 

is  unlikely. 

116 

In  Cd,  the  spin  of  the  state  at  1644  keV,  previous¬ 
ly  given  a  tentative  assignment  of  (2,3)  has  been  shown  to 
be  2+.  A  new  spin  assignment  has  been  shown  to  be  con¬ 
sistent  with  an  assignment  of  3+. 

In  ^^Nd,  the  states  at  1377  and  1471  have  been 
assigned  spin-parities  of  1±  and  2+  respectively.  In  pre¬ 
vious  work,  these  assignments  were  tentative  only.  A  new 
level  at  1787  keV  has  been  proposed. 
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THE  HARMONIC  VIBRATIONAL  MODEL 
FOR  EVEN-EVEN  NUCLEUS 
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Many  models  have  been  formulated  to  explain  the  en- 
ergy  level  schemes  of  nuclei.  The  hydro-dynamic,  or 
liquid  drop,  model  is  a  collective  formalism  that  visual¬ 
izes  the  dynamics  of  the  nucleons  in  the  nucleus  as  mole¬ 
cules  undergiong  motion  in  a  liquid  drop.  In  a  liquid 
drop,  the  lowest  possible  energy  is  achieved  when  the 
drop  has  spherical  symmetry.  An  excited  state  is  said  to 
exist  when  the  surface  of  the  drop  undergoes  deformation. 

The  deformation  of  the  surface  can  be  expressed  as 

R  ( 6  /  $  >  t )  =  R  [1  +  l  a  (t)  Y^(0,<|>)];  |  u  |  <  X  (A-l) 

o  Ay  Ay  X  11 

where  R0  is  the  undeformed  radius  and  X  and  y  refer  to 
the  mode  of  vibration.  The  lowest  mode  of  surface  deforma¬ 
tion  corresponds  to  the  X  =  2  mode.  For  a  nucleus,  the 

largest  mode  of  suface  deformation  corresponds  to 
1/3 

X  ~  A  .  For  higher  modes,  the  vibrations  of  the  nuclear 
surface  are  less  than  the  internucleon  separation,  which 
implies  a  discontinuity  of  the  nuclear  surface.  As  this 
is  unrealistic,  all  higher  modes  of  vibration  are  excluded. 

The  total  deformation  energy  of  a  harmonic  charged 
liquid  drop  is  (Ei70) 

H  -  —  ^  R  I  a  I  ^  —  £  q  I  a  I  2 

H  "  2  XyBX 1  Xy I  +  2  Ay  A  1 aXy 1 

where  B,  represents  the  moment  of  inertia  of  the  drop 
A 


(A-2 ) 


' 
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with  respect  to  changes  in  deformation  and  C-^  is  the 
def ormability  coefficient  which  measures  the  resistance 
of  the  drop  to  deformation.  Introducing  the  variable 

canonically  conjugate  to  ct. 

Ay 

tt,  -  B  a. 

Ay  A  Ay 


gives 


H 


1  £  r  1  i  ,2 

-y  ,  L  —  tta  +  C, 

2  Ay  1  Ay  1  A 


(A- 3) 


This  is  in  the  same  form  as  the  Hamiltonian  for  an  harmonic 
oscillator.  The  solution  is 


E 


(  n^ 
y  Ay 


E  hco,N.  +  -  E 


A 


A  A 


2  A 


hoo  ^  ( 2  A  +  1 ) 


(A-4) 


where  n^  is  the  number  of  oscillators  in  the  Ay  mode, 

N,  is  the  number  of  oscillators  in  the  A  mode; 

A 

N,  =  E  n^  and 

A  y  Ay 

co ^  is  the  characteristic  frequency  of  vibration 
for  the  A  mode.  The  last  term  in  Eq.  (A-4)  is  taken  as 
the  zero  point  energy.  The  characteristic  frequency  co^ 
is  given  by 

(A-5) 

(BI52) 

where  y  is  the  ratio  between  the  Coulomb  and  surface 

2 

tension  energies.  It  is  given  by  y  ~  0.0474  Z  /A. 
a  is  the  surface  tension  coefficient, 

M  is  the  mass  of  a  single  nucleon, 
and  A  is  the  atomic  number. 


oo 


A 


=  C 


A/B 


=  {A ( A  —  1 )  [ (A+2) 


A 


lOy]  4Trg}^//^ 

2 A+l  3MA 


] 
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By  neglecting  the  Coulomb  energy  term,  one  can  obtain 

a  rough  estimate  of  the  ratios  of  the  characteristic 

frequencies.  For  example,  the  ratio  of  the  frequency 

of  the  octupole  (A  =  3)  mode  to  that  of  the  quadrupole 

(A  =  2)  mode  is  approximately  2.  As  w,  .  ~  2.4,  the 

4/u3 

A  -  4  mode  can  be  ignored  when  only  the  low-lying  states 
are  of  interest. 

A  Ay  mode  phonon  is  a  boson  with  a  spin  of  A,  a 
z-component  of  angular  momentum  of  yh,  and  an  intrinsic 
parity  of  (-1)\  When  there  are  two  or  more  phonons 
present,  the  individual  angular  momenta  can  be  coupled 
to  form  eigenstates  of  J  and  J  For  example,  two 
A  =  2  phonons,  corresponding  to  the  N  =  2  state,  form 
a  triplet  of  states  with  spin-parities  of  0+,  2+,  and 
4+.  Considering  only  the  quadrupole  and  octupole  modes 
of  vibration,  a  highly  degenerate,  equally  spaced  level 
structure  is  predicted  (see  fig.  42).  Small  corrections 
to  the  energy  levels  are  thought  to  be  due  to  the  fraction 
of  the  nuclear  matter  involved  in  the  surface  deforma¬ 
tions  undergoing  rotation. 

Since  this  formulism  deals  with  a  spherical  liquid 
drop,  one  would  expect  it  to  be  valid  only  for  nuclei 
that  have  a  spherically  symmetric  ground  state.  These 
nuclei  are  usually  found  where  N  or  Z  is  close  to  a 
magic  number.  Because  of  this,  the  level  scheme  should 
exhibit  marked  shell  effects.  That  is,  there  should  be 
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Figure  42. 


The  level  scheme  predicted  by  the  harmonic  vibra 
tional  model.  The  degeneracy  of  the  levels  has 
been  removed  to  facilitate  identification.  The 
states  are  denoted  by 

I  X ,  N ,  JM  >  where 

X  =  mode  of  vibration 
N  =  the  number  of  phonons 
J  =  the  total  angular  momentum 
M  =  it's  Z  projection 
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a  mingling  of  the  collective  excited  states  of  the 
vibrational  model  with  the  intrinsic  shell  model  states. 

In  the  liquid  drop  model,  the  parameters  B  and  C, 
are  evaluated  classically  by  assuming  an  irrotational 
and  incompressible  drop.  The  particles  in  the  nucleus, 
however  do  not  form  such  a  homogeneous  system  as  the 
particles  in  closed  shells  behave  differently  than  those 
in  unfilled  shells.  The  vibrational  energy  of  the  nucleus 
is  mainly  due  to  the  motion  of  the  nucleons  in  the  unfilled 
shells.  This  is  readily  evident  when  the  experimental 
values  of  B  and  C9  are  obtained  for  a  large  assortment 
of  nuclei.  The  values  can  be  obtained  from  a  considera¬ 
tion  of  the  energy  and  the  reduced  E2  transition  prob¬ 
ability  B  (E2;  0  -*  2)  for  the  excitation  of  the  first 
excited  state  (Ei70). 

“2  =  «VB2)1/2 

and  B(E2;  0  -  2)  =  OZeR^)2  1/2  <A“6> 

2(B2  2 

The  values  obtained  for  the  inertial  parameters  B2 
obtained  are  consistently  higher  than  the  values  predicted 
by  irrotational  flow.  While  in  the  region  of  closed 
shells  the  surface  tension  parameter  C2  is  much  larger 
than  the  liquid  drop  value,  the  C2  parameters  for  nuclei 
with  many  particles  in  unfilled  shells  are  much  lower 
than  the  predicted  values.  This  corresponds  to  the  approach 
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to  instability  of  the  spherical  shape. 

When  the  instability  of  the  nucleus  is  such  that  it 
takes  only  a  small  amount  of  energy  to  break  apart  the 
nucleus,  induced  fission  is  expected.  The  vibrational 
model  gives  the  limit  of  stability  against  fission.  In 
equation  A-5,  it  can  easily  be  seen  that  the  frequency 
for  vibration  because  imaginary  when  the  ratio  between 
the  Coulomb  energy  and  the  surface  energy  approaches  a 
critical  value  y  .  When  this  occurs,  the  negative 
Coulomb  energy  exceeds  the  surface  energy  and  the  nucleus 
is  no  longer  stable  against  surface  deformation.  As 
quadrupole  is  the  lowest  mode  of  vibration,  the  limit 
for  stability  predicted  by  the  vibrational  model  is 


Z2/a  =42.2 


(A- 7) 


Nuclei  approaching  this  limit  have  a  small  deformation 
frequency  to  Therefore  the  amplitude  of  oscillation 
becomes  large,  resulting  in  the  breaking  up  of  the  nucleus. 

Other  fundamental  characteristics  of  vibrational 
excitations  are  an  enhanced  E2  transition  from  the  first 
excited  state  to  the  ground  state  as  compared  to  single 
particle  estimates  and  a  small  static  electric  quadrupole 
moment.  The  vibrational  model  also  imposes  certain 
restrictions  on  the  type  of  y  decay  that  can  occur.  As 
the  second  2+  excited  state  results  from  the  coupling  of 
two  X  =  2  phonons,  the  decay  of  this  level  to  the  ground 


state  is  forbidden.  The  decay  can  proceed  only  via  a 
cascade  to  the  first  excited  state.  In  this  cascade. 
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Ml  radiation,  which  is  preferred  by  single  particle 
motion,  is  strictly  forbidden;  the  transition  taking 
place  only  through  E2  radiation. 


